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Abstract 

A complete classification of the perfect binary one-error-correcting codes of length 15 as well as their extensions 
of length 16 is presented. There are 5 983 such inequivalent perfect codes and 2 165 extended perfect codes. Efficient 
generation of these codes relies on the recent classification of Steiner quadruple systems of order 16. Utilizing a 
' result of Blackmore, the optimal binary one-error-correcting codes of length 14 and the (15, 1024,4) codes are 

also classified; there are 38 408 and 5 983 such codes, respectively. 
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I. Introduction 

Consider the space FI> 1 of dimension n over the Galois field F 2 = {0, 1}. A binary code of length 
q ; n is a subset of W%. The (Hamming) distance d(x, y) between two codewords x, y is the number of 
1— '> coordinates in which they differ, and the (Hamming) weight wt(x) is the number of nonzero coordinates. 
(*C) ! The support of a codeword is the set of nonzero coordinates, supp(x) = {i : X{ ^ 0}. Accordingly, 
d(x,y) = wt(x-y) = |supp(x-y)|. 

A code has minimum distance d if d is the largest integer such that the distance between any dis- 
lO ■ tinct codewords is at least d. Then the balls of radius [(d — l)/2j centered around the codewords are 
! nonintersecting, and the code is said to be a \_(d — 1)/2J -error-correcting code. If these balls tile the 
^■O whole space, then the code is called perfect. The parameters of perfect codes over an alphabet of prime 
■ order are well known [1], and perfect binary codes exist with d — 1; d — n; d — (n — l)/2 for odd n; 
O ■ d = 3,n = 2 m — 1 for m > 2; and d — 7, n — 23. The first three types of codes are called trivial, the 
j> ! fourth has the parameters of Hamming codes, and the last one is the binary Golay code. A perfect code 
k> ' with minimum distance d is also called a \_(d — 1)/2J -perfect code. 

A binary code with length n, minimum distance d, and M codewords is called a (n, M, d) code. In this 
. 5^ ! notation a binary 1-perfect code is a (2 m — 1, 2 2m ~ m ~ 1 , 3) code. Two related families are the extended and 
shortened 1-perfect codes, which have parameters (2 m , 2 2m_m_1 , 4) and (2 m — 2, 2 2m ~ m ~ 2 , 3), respectively. 

Existence of binary 1-perfect codes follows from the existence of Hamming codes, which are the unique 
linear 1-perfect codes. Still constructing all 1-perfect codes is a longstanding open problem. It makes sense 
to approach this issue by considering the number of inequivalent codes (or more formally the number of 
equivalence classes). Two codes are said to be equivalent if one is obtained from the other by permuting 
coordinates and adding a constant vector; a formal definition appears in Section lU 

There is trivially a unique 1-perfect code of length 3. Zaremba [2] showed that also the 1-perfect code 
of length 7 is unique. However, already the next case of length 15 has until this work withstood all 
attempts of complete classification, although several constructions of such codes have been published; see 
the surveys [3], [4]. It turns out that these results were not far from a complete classification as for the 
number of codes found. The growth of the number of 1-perfect binary codes is double exponential in the 
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length of the code, see [5] for a lower bound on this number. For an in-depth treatment of the topic of 
classifying combinatorial objects, see [6]. 

The aim of the current work is to obtain a complete classification of inequivalent 1 -perfect binary 
codes of length 15. By computer search it is here shown that their number is 5 983. Also the codes 
obtained by extending, shortening or extending and shortening are classified; the numbers of (16, 2 048, 4), 
(14, 1 024, 3) and (15, 1 024, 4) codes turn out to be 2 165, 38 408 and 5 983 respectively. 

In the rest of the paper we document the classification of the extended 1 -perfect codes of length 16, 
which yields classifications of the 1 -perfect codes of length 15 and the shortened 1 -perfect codes of 
length 14. In Section [XT] we define some concepts and consider construction of extended 1 -perfect codes 
via Steiner systems. In Section [III] we present algorithms for detecting and rejecting equivalent codes, 
and in Section [TV] we take a brief look at the results; a separate, more detailed study of the classified 
codes will appear in a separate paper [7]. Finally, in Section |V] we give a consistency check for gaining 
confidence in the computational results. 

II. Preliminaries and Construction 

A permutation n of the set {1,2,..., n} acts on codewords by permuting the coordinates: 7r((ci, C2, . . . , c n )) = 
(Ctt-ih), 0^-1(2), • • • , Cjr-ifa)). Pairs (n, x) form the wreath product S2lS n , which acts on codes as (it, x)(C) = 
7r(C + x) = 7r(C) +7r(x). Two codes, C\ and Ci, are isomorphic if C\ = n(C2) for some n and equivalent 
if C\ = 7r(C 2 + x) for some n, x. 

The automorphism group of a code C, Aut(C), is the group of all pairs (n, x) such that C = 7r(C + x). 
Two important subgroups of Aut(C) are the group of symmetries, 

Sym(C7) = {tt : n(C) = C} 

and the kernel 

Ker(C7) = {x : C + x = C}. 

If the code contains the all-zero word, 0, then the elements of the kernel are codewords. 

A Steiner system S(t, k, v) can be viewed as a code S C F2 with the property that each codeword of S 
has weight k, and for any y G with wt(y) = t, there is a unique x G S such that supp(y) C supp(x). 
Usually Steiner systems are defined as set systems rather than codes, but our definition is more directly 
applicable for this work. The parameter v is the order of the system. Steiner systems 5(2, 3, v) and 
5(3, 4, v), which are called Steiner triple systems and Steiner quadruple systems, respectively, are related 
to 1 -perfect codes in the following way. If C is a 1 -perfect binary code of length v and x G C, then the 
codewords of C + x with weight 3 form a Steiner triple system of order v. Similarly, if C is an extended 
1 -perfect binary code and x G C, then the codewords of C + x with weight 4 form a Steiner quadruple 
system. These systems are, respectively, the neighborhood triple system and neighborhood quadruple 
system associated with the code and the codeword. 

As a starting point for the classification of the extended 1 -perfect binary codes of length 16, we have 
the classification [8] of Steiner quadruple systems of order 16; there are 1 054 163 such designs. We want 
to find, for each 5(3,4,16), all extended 1-perfect binary codes in which it occurs. This can be done 
by puncturing any coordinate, augmenting the resulting code to 1-perfect codes in all possible ways, and 
finally extending every resulting code with a parity bit. 

When augmenting a set of codewords to a 1-perfect code, we consider a 1-perfect code as a set of balls 
with radius one that form a partition of the ambient space. Accordingly, finding a code (with specified 
codewords) is a special case of the exact cover problem, where we are given a set 5 and a collection 
U of its subsets, and the task is to form a partition of 5 by using sets in U. Let the set Q contain the 
codewords obtained by puncturing the all-zero codeword and its neighborhood quadruple system. In this 
case we have S = ¥ 2 5 \ B(Q), and U = {£(x) : £(x) n B(Q) = 0}, where 5(x) = {y : d(x, y) < 1} 
and B(C) = {B(x) : x G C}. We use the libexact software [9] for solving such instances of the exact 
cover problems. In the search we could in fact have made use of the fact that all 1-perfect binary codes 
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are self-complementary — in other words, the all-one word is always in the kernel — but this would not 
have had any practical significance as the search was rather fast. 

Since the all-zero word and its neighborhood quadruple system contain 141 of the 2 048 codewords, 
1907 new codewords are needed. Searching for these was a remarkably easy computational task; on 
average the search trees in which codes were found had 1978 nodes and those in which no codes were 
found had 3 nodes. 

III. ISOMORPH REJECTION 

The general framework by McKay [10] was used to carry out isomorph rejection, although a less 
sophisticated method would have sufficed in this work. 

Recall that we augment a Steiner quadruple system Q to an extended 1 -perfect code C. We accept C 
if it passes the following two tests; otherwise it is rejected. First we require that C shall be the minimum 
(with respect to some practically computable total order of codes) under the action of Aut(Q). Second, we 
compute the canonical equivalence class representative c E (C), consider ir, x for which 7r(C + x) = ce(C) 
and require that x and are on the same Aut(C) orbit (we define ce so that xgC always holds). 

When the extended 1 -perfect codes have been classified, classifying the 1 -perfect codes is straightfor- 
ward. All 1-perfect codes are obtained by puncturing the extended codes, and the resulting 1-perfect codes 
are equivalent if and only if they are obtained by puncturing the same extended code at coordinates which 
are in the same orbit of the automorphism group. 

A complete classification of the (14,1024,3) codes is obtained similarly, since each such code is 
obtained by shortening a unique (up to equivalence) 1-perfect code of length 15; this result was proved 
by Blackmore [11]. Although a code can be shortened at any coordinate in two ways, by selecting the 
codewords with or 1 in a certain coordinate, both selections lead to equivalent codes. This follows from 
the fact that every 1-perfect binary code is self-complementary. 

Furthermore we note that any (15, 1 024, 4) code is obtained by extending a (14, 1 024, 3) code with a 
parity bit. Hence all such codes are obtained by shortening and extending a perfect code, or equivalently 
removing all words of chosen parity. As the perfect codes are self-complementary, we get (up to equiva- 
lence) same code by chosing either odd or even parity. As this mapping is reversible, we conclude that 
there is one-to-one correspondence between equivalence classes of (15,2 048,3) codes and equivalence 
classes of (15, 1 024, 4) codes, and in both cases their number is 5 983. 

Let C be a (15, 1024,4) code C and let C be the corresponding perfect code. The group Aut(C') 
contains Aut(C) as a subgroup, and Aut(C') has one more generator than Aut(C), namely the all-one 
codeword. Accordingly |Aut(C")| = 2|Aut(C)|. 

We still have to describe an algorithm for canonical labeling. The most straightforward approach of 
using the general purpose isomorphism nauty [12] is rather slow on codes as large and regular as the 
(16,2 048,4) codes; this was also noted by Phelps [13]. Hence a tailored approach is necessary. The 
method presented below has a lot in common with the one described in [13]. An alternative method based 
on minimum distance graphs would also work [14], cf. [15]. 

A triangle consists of 3 codewords with mutual distance 4. Triangles constitute an easily computable 
and rather sensitive invariant of Steiner quadruple systems. Distinguishing the isomorphism classes of the 
neighborhood quadruple systems of a code also constitutes an invariant of the extended 1-perfect codes. 
These two invariants turned out to be useful for speeding up our computations. 

A canonical isomorphism class representative C/(C) for a code C can be computed by using nauty 
to label a corresponding graph canonically. Moreover, nauty computes generators of the group Sym(C). 
Canonical equivalence class representative can be defined as c E (C) = min{c/(C + x) : x e C}, where 
the minimum is again taken with respect to some practically computable total order of codes. 

Note that two codewords, x and y, are in the same orbit of Aut(C) if and only if cj(C+x) = c/(C+y). 
Because of this, if we know that x and y are in the same orbit, and cj(C + x) has been computed, then 
there is no need to compute cj(C + y). Also if cj(C + x) = cj(C), then nauty yields a permutation ir 
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such that 7r(C + x) = C. The pairs (-7T, x) are coset representatives of Aut(C) with respect to Sym(C), 
so we get generators of the group Aut(C). 

IV. Results 

There are exactly 2 165 inequivalent extended 1 -perfect codes of length 16, 5 983 inequivalent 1 -perfect 
codes of length 15, 38 408 shortened 1-perfect codes of length 14 and 5 983 (15, 1 024, 4) codes. The orders 
of the automorphism groups of the codes are presented in Tables HI [H] and HHl As noted in Section Hill the 
order of the automorphism group of a (15, 1 024, 4) code is half of the order of the automorphism group 
of the corresponding perfect code. 

The codes have been made available in electronic form by including them in the source of the arXiv ver- 
sion of this paper. Downloading |http : //arxiv. org/e-print /0806 . 2513v3| and uncompressing 
it with gunzip and tar yields the files perfect 15 and extendedl 6. 

Only 15 590 of the 1 054 163 nonisomorphic S(3, 4, 16) can be augmented to a 1-perfect code, and the 
total number of extensions is 22 814. The computationally intensive part of this result was the earlier 
classification of £(3,4, 16), which required several years of CPU time, while all searches described in 
this paper took only a couple of hours of CPU time. 

A detailed study of the properties of the classified codes will appear in a second part of this article [7]. 

TABLE I 



Automorphism groups of (16,2 048,4) codes 



|Aut(C)| 


# 


|Aut(C)| 


# 


|Aut(C)| 


# 


128 


11 


5 376 


1 


196 608 


6 


192 


5 


6144 


23 


262 144 


3 


256 


105 


8192 


174 


344 064 


1 


384 


9 


10752 


2 


393 216 


3 


512 


377 


12288 


22 


524 288 


2 


672 


2 


16384 


103 


688 128 


1 


768 


19 


24576 


12 


786432 


2 


1024 


416 


32 768 


47 


1 572 864 


3 


1344 


1 


43 008 


2 


2359 296 


1 


1536 


21 


49152 


18 


2752512 


1 


1920 


1 


61440 


1 


3 145 728 


1 


2048 


394 


65 536 


33 


5 505 024 


2 


2 688 


1 


86016 


3 


6291456 


1 


3 072 


18 


98 304 


12 


660 602 880 


1 


4096 


298 


131072 


6 







V. Consistency check 

To get confidence in the results, we performed a consistency check similar to the one used, for example, 
in [8]. In this check we count the total number of codes in two different ways and ensure that the results 
agree. 

First we consider the set C of equivalence class representatives obtained in the classification. By the 
orbit- stabilizer theorem, the total number of extended 1-perfect codes is 

16! ■ 2 16 



^Aut(C)' 

where 16! • 2 16 is the order of the wreath product group acting on the codes. 

Let Q consist of the representative Steiner quadruple systems, and let E(Q) be the number of all 
extended 1-perfect codes obtained by augmenting Q. Applying the orbit-stabilizer theorem, we get the 
expression 

1 ^ 16! • 2 16 • E{Q) 



E 



2 048^ Aut(Q) 
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TABLE II 

Automorphism groups of (15, 2 048, 3) codes 



|Aut(C)| 


# 


|Aut(C)| 


# 


|Aut(C)| # 


8 


3 


512 


1017 


24576 7 


12 


3 


672 


3 


32 768 8 


16 


5 


768 


32 


43 008 4 


24 


10 


1024 


697 


49 152 10 


32 


138 


1536 


17 


65 536 5 


42 


2 


2048 


406 


98 304 1 


48 


12 


2688 


1 


131072 1 


64 


542 


3 072 


37 


172032 1 


96 


22 


3 840 


1 


196 608 5 


120 


1 


4096 


202 


344064 2 


128 


1230 


5 376 


4 


393 216 2 


192 


18 


6144 


35 


589 824 1 


256 


1319 


8192 


94 


41287 680 1 


336 


3 


12 288 


7 




384 


30 


16384 


44 





TABLE III 

Automorphism groups of (14, 1 024, 3) codes 



|Aut(C)| 


# 


|Aut(C)| 


# 


|Aut(C)| 


# 


1 


5 


168 


1 


8192 


80 


2 


75 


192 


80 


12 288 


18 


3 


8 


256 


4 392 


16384 


14 


4 


425 


336 


5 


21504 


1 


6 


39 


384 


114 


24576 


15 


8 


1 162 


512 


2469 


32768 


14 


12 


56 


768 


30 


49152 


1 


16 


3 465 


1024 


1346 


65 536 


1 


21 


4 


1344 


1 


86016 


1 


24 


39 


1536 


54 


98 304 


2 


32 


7311 


2048 


527 


172 032 


1 


48 


59 


2 688 


6 


196 608 


2 


64 


9 068 


3 072 


55 


1376 256 


1 


96 


49 


4096 


222 






128 


7 172 


6144 


18 







where the division by 2 048 is necessary since each code is counted once for each codeword. Both 
formulas yield the same result, 2 795 493 027033 907200. Similarly we also counted the 1-perfect codes 
and shortened 1-perfect codes in two different ways; their number is 1397 746 513 516 953 600. Indeed, 
there are twice as many extended 1-perfect codes as there are 1-perfect codes, since each 1-perfect code 
admits two extensions: one with even parity bit and one with odd. Similarly, we get a bijection from the 
1-perfect codes to shortened 1-perfect codes if we shorten each code by taking, for instance, the codewords 
with value in coordinate 15 and removing that coordinate. Thus there are equally many 1-perfect codes 
and shortened 1-perfect codes. 
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